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A Lagrangian density is provided, that allows to recover the Z4 evolution system from an action 
principle. The resulting system is then strongly hyperbolic when supplemented by gauge conditions 
like T+log' or 'freezing shift', suitable for numerical evolution. The physical constraint = can 
be imposed just on the initial data. The corresponding canonical equations are also provided. This 
opens the door to analogous results for other numerical-relativity formalisms, like BSSN, that can 
be derived from Z4 by a symmetry-breaking procedure. The harmonic formulation can be easily 
recovered by a slight modification of the procedure. This provides a mechanism for deriving both 
the field evolution equations and the gauge conditions from the action principle, with a view on 
using symplectic integrators for a constraint-preserving numerical evolution. The gauge sources 
corresponding to the 'puncture gauge' conditions are identified in this context. 
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I. INTRODUCTION 



The role of action principles is so crucial in theoretical 
physics that its importance can not be overemphasized. 
In the case of the General Relativity, the standard ac- 
tion was proposed by Hilbert since the very beginning 
of the theory, although the Hamiltonian formulation had 
to wait for decades [J, [|J . The reason for this long de- 
lay is probably related to the complexity of the Cauchy 
problem for Einstein's equations, which becomes mani- 
fest in the 3+1 (space plus time) decomposition Q. The 
coordinate gauge freedom produces a mismatch between 
the number of dynamical fields and that of true evolution 
equations: four of the field equations are indeed (energy- 
momentum) constraints. This rich structure opens the 
door to many different approaches. 

On the other hand, by the end of the past century, 
some hyperbolic extensions of Einstein's equations were 
developed with a view on numerical relativity applica- 
tions [3-6]. This emergent field is now more mature: 
there are two main formalisms currently used in nu- 
merical simulations. One is BSSN 0, Q, working at 
the 3+1 level, and the other is the class of general- 
ized harmonic formalisms [ij-fllj]. working at the four- 
dimensional level. A unifying framework is provided by 
the Z4 formalism pjj], which allows to recover the gen- 
eralized harmonic one by relating the additional vector 
field Z^ with the harmonic 'gauge sources' @. On the 
other hand, it allows to recover (a specific version of) 
BSSN by a symmetry-breaking process in the transition 
from the four-dimensional to the three-dimensional for- 
mulations [HEi]. 

There is a growing interest in incorporating 
the new hyperbolic formulations into the La- 
grangian/Hamiltonian framework. An example is 
the usage of the 'densitized lapse' Q as a canonical 
variable, leading to a modification in the standard form 
of the canonical evolution equations 15]. Reciprocally, 



there are very recent attempts of modifying the ADM 
action in order to incorporate coordinate conditions 
of the type used in numerical relativity [lril \vf\ . with 
a view on using symplectic integrators for the time 
evolution, which could ensure constraint preservation in 
numerical simulations [l8j]. On a different context, a well 
posed evolution formalism developed from a Lagrangian 
formulation could be a good starting point for Quantum 
Gravity applications. 

In this paper we derive for the first time the Z4 for- 
malism from an action principle by introducing a La- 
grangian density which generalizes the Einstein-Hilbert 
one. We also provide the corresponding Hamiltonian, 
via the Legendre transformation. This is a crucial step 
towards the Hamiltonian formulation of other numerical- 
relativity formalisms, like BSSN. On the other hand, we 
recover the generalized harmonic formulations as usual, 
by relating the additional vector field Z M with the har- 
monic 'gauge sources'. This mechanism is generalized, 
by identifying the gauge sources which correspond to the 
current numerical-relativity coordinate conditions, as we 
show explicitly for the 'puncture gauge': the combination 
of the ' 1+log' slicing and the gamma-driver conditions. 



II. THE ACTION PRINCIPLE 



Let us consider the generic action 



S = 



C 



(1) 



with a Lagrangian density which generalizes the Einstein- 
Hilbert one by including an extra four- vector Z^, namely 

C = ^/g g^[Ruv + 2VuZ v ] (2) 

(we restrict ourselves to the vacuum case) , with the Ricci 
tensor written in terms of the connection coefficients 



R. 
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(round brackets denote symmetrization). 

Now let us follow the well-known Palatini approach, by 
considering independent variations of the metric density 
= \/9 9^" i the connection coefficients V and the 
vector Z^. From the W variations we get directly the 
Z4 field equations [l2| 



(4) 



which are currently used in many numerical-relativity de- 
velopments [l4| . 

From the and the variations we get a coupled 
set of equations which can be decomposed in a covariant 
way into the tensor equation 



v p <r = 0. 



(5) 



which fixes the connection coefficients in terms of the 
metric, and the vector condition 



Z IM = 0. 



(6) 



Let us note here the different role of the conditions (0 
and ([6]) . As there are much more independent connection 
coefficients than evolution equations in we will con- 
sider condition ([5]) as a constraint enforcing the metric 
connection 'a posteriori', that is after the variation pro- 
cess. In this way, we will ensure that equations Q are 
identical to the original Z4 equations, rather than some 
affine generalization. For this reason, we will assume a 
metric connection everywhere in what follows. 

The case of condition (jSJ) is different, as the Z4 equa- 
tions ((3]) actually provide evolution equations for every 
component of Z^. Then, ([6]) is a standard primary con- 
straint and we have a choice among different strategics for 
dealing with it. If we enforce ^ into the Z4 field equa- 
tions (HI, we get nothing but Einstein's equations. This is 
not surprising because our Lagrangian obviously reduces 
to the Einstein-Hilbert one when Z^ vanishes. The prob- 
lem is that the plain Einstein field equations do not lead 
directly to a well-posed initial data problem. This is why 
the original harmonic formulation [2H - [23t was used in- 



stead in the context of the Cauchy problem [24| . For the 
same reason, other formulations (BSSN @,Hj], generalized 
harmonic [9Hli"j|. Z4 [H|,[l3[) are currently considered in 
numerical relativity. 

We can alternatively follow a different strategy. In- 
stead of enforcing ^ , we can deal with this condition as 
an algebraic restriction to be imposed just on the initial 
data, that is 



Za \t=o = 



(7) 



In this way, we can keep the Z4 field equations which, 
when supplemented with suitable coordinate conditions, 
lead to a strongly hyperbolic evolution system [IH, [l4| ■ 
The consistency of this 'relaxed' approach requires that 
the constraint ^ should be actually preserved by the Z4 
field equations ((U). In this way, the solutions obtained 



from initial data verifying (J6j) will actually minimize the 
proposed action (JTJ) . 

Allowing for the conservation of the Einstein tensor, 
which is granted after the metric connection enforcement, 
we derive from Q the second-order equation, linear- 
homogeneous in Z 



V„ [ V»Z U + V V Z» - (V p Z p ) g" v } = . 



(8) 



It follows that the necessary and sufficient condition for 
the preservation of the constraint Q is to impose also its 
first time-derivative conditions in the initial data, that is 



(do 



Z„)|t=o = 0. 



(9) 



Note that, allowing for ([7]) and the Z4 field equations, 
the secondary constraints ^ amount to the standard 
energy and momentum constraints, which are then to be 
imposed on the initial data in addition to ([7]). 

This ' relaxed' treatment of the constraints (JBJ) may 
look unnatural. But it is just the reflection of a com- 
mon practice numerical relativity ('free evolution' ap- 
proach) , where four of the ten field equations (the energy- 
momentum constraints) arc not enforced during the evo- 
lution, being imposed just in the initial data instead. The 
introduction of the extra four-vector in the Z4 formalism 
actually provides a simpler implementation of the same 
idea. 



III. HAMILTONIAN FORMALISM 

A detailed look at the Lagrangian density (|2|) shows 
that the time derivatives of most of the variables are not 
present in C. The only exceptions are the combinations 

+ (to) 

This dynamical subset of variables can be decomposed 
into 

{t%, z,- 1 - (r^.-r° 0i ), Z -±T k k0 } (ij,k = 1,2,3), 

The corresponding canonical momenta are given, respec- 
tively, by 

{ h ij , 2h 0i , 2h 00 }. (12) 
The remaining quantities 

{rV,r° 0/J ( i3) 

can be considered as a sort of Lagrange multipliers, in- 
troducing constraints into the dynamical system, as their 
time derivatives do not enter the Lagrangian (see for in- 
stance ref. [Hj]). 

Note that our Lagrangian @ is linear in the time 
derivatives. This means that the relationship between 
field velocities and momenta can not be inverted. The 
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canonical momenta corresponding to (fl2l [TBI vanish 
identically, and the ones corresponding to (fTTj) coincide 
with the metric density components (fT2|) . This means 
that the Lagrangian @ is a singular one: the correspond- 
ing (constrained) canonical formalism can be developed 
following the work of Dirac 20] . 

We will rather sketch here a simpler approach, by per- 
forming a limited Legendre transform, in the sense that 
it will only affect the dynamical subset (fTTj) , with canoni- 
cal momenta ()12j) . The remaining quantities will be con- 
sidered as Lagrange multipliers, for which no Legendre 
transformation is required. We obtain in this way the 
Hamiltonian function 



H = -h^{ d k [ T% - (T% 
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(14) 



where the metric densities are considered here as the 
canonical momenta associated to the dynamical fields 

(UJ. 

The Hamilton equations for the fields (fTTj) are precisely 
the Z4 equations ((4]). The Hamilton equations for their 
momenta (1121) can be written as 



dp hi* 



= - V 

= - W r* 



h ij {T p pQ -2Z )- 2h p{ - i Y j) 



pO 



(15) 
(16) 
(17) 



which must be supplemented with the constraints derived 
from the Lagrange multipliers subset (fT5|) . A straight- 
forward calculation shows that the full set of Hamilton 
equations is still equivalent to the Z4 equations (jT]), plus 
the metric connection requirement ([5]), plus the vanish- 
ing of (JSJ). Indeed, allowing for (0 [HJ) , the subsystem 
(|15H17p is verified identically. 

Note that in all our developments we have preserved 
general covariance. Our action integral (UJ) is a true scalar 
and, in spite of other alternatives, we have avoided the 
addition of total divergences which could have simplified 
our developments to some extent, at the price of adding 
boundary terms. This means that we keep at this point 
the full coordinate-gauge freedom. 

This is reassuring from the theoretical point of view, 
but it can be a disadvantage if one is planning to use 
symplectic integrators for numerical evolution, as the re- 
quired coordinate conditions must be supplied from the 
outside of the canonical formalism. This is why some 
recent works are trying to incorporate the coordinate 
conditions, via Lagrange multipliers, into the canonical 
framework [III [171. 



IV. GENERALIZED HARMONIC SYSTEMS 

There is still another possibility, which allows a more 
direct specification of a coordinate gauge at the price of 
breaking the covariance of the evolution equations, which 



are currently used in many numerical-relativity develop- 
ments (l4j ]. We can enforce in the Z4 equations (01 the 
following assignment for 



Z p 



- T p q pa 
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pp 



(18) 



The vanishing of Z^ will amount in this way to the 'har- 
monic coordinates' condition, which can be considered 
then as a constraint to be imposed just in the initial 
data, that is 



pa\ 







(19) 



(note that the extra field Z^ has disappeared in the pro- 
cess). The resulting field equations 



pp p 

x pi/ x p 







(20) 



lead, after imposing the metric connection condition ([5]) , 
to the manifestly hyperbolic second-order system 

9 p °d% 9»u = 2 9 pa 9 aP [d a5 pp d p9al ,~T wa T w p] . (21) 

This corresponds to the classical harmonic formulation 
of General Relativity |2l| - |23| . which is known to have a 
well-posed Cauchy problem [24j . 

We have derived in this way the harmonic formalism 
through the non-covariant prescription (|18[) . The har- 
monic constraint (|T9l is automatically preserved by the 
resulting (harmonic) evolution system, provided that we 
also enforce the energy-momentum constraints on the 
initial data. This can be seen in a transparent way 
by replacing directly (TT81 into the covariant constraint- 
evolution equation (|8j and then into the resulting condi- 
tions ©. 

The prescription (|18p can be generalized in order to 
enforce other coordinate gauges that are also currently 
used in numerical relativity. The simpler formulations 9, 
[TlT | correspond to the replacement 



Z p = -- (T p 
2 V 



£P), 



(22) 



where the ' gauge sources' H p are explicit functions of the 
metric and/or the spacetime coordinates. More general 
choices of H p , like that of ref . [25[ , would require a more 
elaborate treatment. 

The same mechanism can be applied to coordinate con- 
ditions derived in the 3+1 framework, where the space- 
time line element is decomposed as 



ds 1 



dt 2 + jij (dx l + F dt) (dx j + /3 j dt) . (23) 



The spacetime slicing is given by the choice of the time 
coordinate. In this context, the harmonic slicing condi- 
tion can be generalized to [l3| 



(d t -l3 k d k ) a = -fa 2 trK, (24) 

where K%j — -aT^ stands for the extrinsic curvature 
of the time slices. The case / = 1 corresponds to the 
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r nnn = l/a 2 (d t - I3 r dr) a T nnk = -dklna 

V knn = l/a 2 j kj (d t - P r dr)P 3 + d k lna V ni j = -K is 

r m = Ki, - i/a 7lfe dj p k Vkij = (3) r kij 



TABLE I: The 3+1 decomposition of the four-dimensional 
connection coefficients. The index n is a shorthand for the 
contraction with the unit normal n M . 

harmonic time-coordinate choice, whereas the choice / = 
2/q corresponds to the popular ' 1+log' time slicing. 

In order to get the replacement, of the form (f2"2"j) . which 
connects this condition with our formulation, we must 
rewrite (|24l) in a four-dimensional form with the help of 
the unit normal n p to the constant time hypersurfaces, 
that is 

n^ = aS° p = + (25) 

We can now decompose the four-dimensional Christoffel 
symbols in terms of the standard 3+1 quantities (see Ta- 
ble I). This provides a convenient way of translating 3+1 
conditions like (j2"4"]l in terms of four-dimensional objects. 

We can obtain in this way, after an straightforward 
calculation, the gauge sources corresponding to the class 
of slicing conditions (l24l) . namely 

H° = (1 - l/f) T a pa n"n- . (26) 

We will use now (|22[) for replacing the quantity Z° in the 
Z4 equations. Its evolution equation gets transformed in 
this way into a second order evolution equation for the 
lapse function a, which governs the spacetime slicing. As 
the first-order slicing condition ([24]) has been translated 
into a specification of Z°, and allowing for ©, (|24p will 
become a first integral of the second order evolution sys- 
tem: we can impose it just in the initial data together 
with the energy-momentum constraints. This approach 
is new in 3+1 formalisms, but a common practice in the 
harmonic-like ones. 

The same technique can be used for ' gamma-driver' 
shift prescriptions. A first-order reduction of the original 
' gamma-freezing' condition [2fj| is given by [27j 

(d t - p k 3 k ) f3 i = n P - V (3\ (27) 

where T l stands here for the contraction of the three- 
dimensional conformal connection, that is 

f i = 7 «V(r J „ + i r rsj ). (28) 

The corresponding ' gamma-driver' gauge sources are 
given by 

a 2 1 
Hi = (1 ) Y lpa nPn° + - Y pal gP a (29) 

1 a 2 

+ (q )Y pal n p n a - j] In g i . 



We can use again (|22|) . this time for replacing the space 
vector Zi in the Z4 equations. Its evolution equation get 
transformed in this way into a second order evolution 
equation for the shift components /3 1 , which determine 
the time lines. Again, the first-order gamma-driver con- 
dition (|27p becomes a first integral of the resulting (sec- 
ond order) shift evolution equation. At the same time, 
one gets rid of the additional variables Zi (as we did for 
Z° with the analogous replacement, leading to the lapse 
evolution equation). 



V. CONCLUSIONS AND OUTLOOK 

In summary, we are proposing the action ([1]), which 
generalizes the Einstein-Hilbert one. Starting from this 
action one gets directly the Z4 field equations, plus the 
metric connection condition (which is to be enforced ' a 
posteriori' in our Palatini approach), plus the constraints 
© stating the vanishing of Z p . We have shown how 
a suitable treatment of these constraints allows working 
with the Z4 covariant evolution in the way one usually 
does in numerical relativity. The price to pay for this 
general-covariant approach is that closing the evolution 
system requires a separate coordinate gauge specification. 
The challenge is then to incorporate the evolution equa- 
tions for the gauge-related quantities (lapse and shift) 
into the canonical formalism, either via Lagrange multi- 
pliers [ljl [ItJ or by any other means. 

We have also presented an alternative strategy, based 
in the ' gauge sources' approach, which characterizes the 
generalized harmonic formalisms. This allows to dis- 
pose of the additional Z p vector field by enforcing at 
the same time the required coordinate conditions by 
means of some generalized gauge sources. The advan- 
tage of this second approach, at the price of getting a 
non-covariant evolution system, is that it can allow a 
direct use of symplectic integrators, devised to ensure 
constraint preservation during numerical evolution (see 
for instance ref. [Hj]). We have actually identified the 
gauge sources corresponding to some standard 3+1 gauge 
conditions, like the 'puncture gauge' consisting of the 
' 1+log' lapse plus the gamma-driver shift prescriptions. 
The fact that these popular gauge conditions can fit into 
a Lagrangian/Hamiltonian approach, in the way we have 
shown, opens the door to new numerical relativity devel- 
opments. 

Acknowledgments The authors are indebted with 
Luis Lehner for useful comments and suggestions. This 
work has been jointly supported by European Union 
FEDER funds and the Spanish Ministry of Science 
and Education (projects FPA2007-60220 and CSD2007- 
00042). C. Bona-Casas acknowledges the support of the 
Spanish Ministry of Science under the FPU program. 
C. Palenzuela acknowledges support from Louisiana 
State University through the NSF grant PHY-0803629. 



5 



[1] P.A.M. Dirac, Proc. Roy. Soc. A246, 333 (1958). 

[2] R. Arnowitt, S. Deser and C.W. Misner, in Gravitation: 

an introduction to current research, ed. L. Witten (Wiley: 

New York 1962). |gr-qc/0405109| 
[3] Y. Choquet-Bruhat and T. Ruggeri, Comm. Math. Phys. 

89, 269 (1983). 
[4] C. Bona and J. Masso, Phys. Rev. Lett. 68 1097 (1992) 
[5] C. Bona, J. Masso, E. Seidel and J. Stela, 

Phys. Rev. Lett. 75 600 (1995). 
[6] A. Abrahams, A. Anderson, Y. Choquet-Bruhat and 

J. W. York, Phys. Rev. Lett. 75, 3377 (1995). 
[7] M. Shibata and T. Nakamura, Phys. Rev. D52 5428 

(1995). 

[8] T. W. Baumgarte and S. L. Shapiro, Phys. Rev. D59 

024007 (1999). 
[9] H. Friedrich, Comm. Math. Phys. 100, 525(1985). 
[10] F. Pretorius, Phys. Rev. Lett. 95 121101 (2005). 
[11] L. Lindblom et al, Class. Quantum Grav. 23, S447 
(2006). 

[12] C. Bona, T. Ledvinka, C. Palenzuela and M. Zacek, 

Phys. Rev. D 67, 104005 (2003). 
[13] C. Bona, T. Ledvinka, C. Palenzuela and M. Zacek, 

Phys. Rev. D 69, 064036 (2004). 
[14] C. Bona, C. Palenzuela and C. Bona-Casas Elements 

of Numerical Relativity and Relativistic Hydrodynamics, 

Lect. Notes Phys. 783 (Springer, Berlin-Heidelberg-New 

York 2009). 

[15] A. Anderson and J.W. York, Phys. Rev. Lett. 81, 1154- 
1157 (1998). 

[16] J. D. Brown in, Quantum Mechanics of Fundamen- 



tal Systems: The Quest for Beauty and Simplicity, 
ed. by M. Henneaux and J. Zanelli (Springer, Berlin- 
Heidelberg-New York 2009). larXiv70 803.0334 [gr-qc]. 
[17] D. Hilditch and R. Richter, larXiv:1002.4TT9V l [gr-qc] 
(2010). 

[18] C. Di Bartolo, R. Gambini and Jorge Pullin, 

J. Math. Phys. 46, 032501 (2005). 
[19] C. Lanczos, The variational principles of Mechanics 

Dover, New York (1986). 
[20] P.A.M. Dirac, Lectures on Quantum Mechanics (Yeshiva 

University Press, New York 1964). 
[21] T. De Donder, La Gravifique Einstenienne Gauthier- 

Villars, Paris (1921). 

The Mathematical Theory of Relativity, (Massachusetts 

Institute of Technology, 1927). 
[22] K. Lanczos, Ann. Phys. 13, 621 (1922). 

Z. Phys. 23, 537 (1923). 
[23] Fock, V.A., The Theory of Space, Time and Gravitation, 

Pergamon, London (1959). 
[24] Y. Choquet (Foures)-Bruhat, Acta Matematica 88, 141 

(1952). 

'Cauchy problem' in Gravitation: An introduction to 
Current Research, ed. by L. Witten, Wiley, (New York 
1962). 

[25] F. Pretorius, Class. Quant. Grav. 23, S529 (2006). 
[26] M. Alcubierre et al, Phys. Rev. D67, 084023 (2003). 
[27] J. R. Van Meter, J. G. Baker, M. Koppitz and D-I. Choi, 
Phys. Rev. D73, 124001 (2006). 



